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Abstract — In cognitive radio systems, one of the main re- 
quirements is to detect the presence of the primary users' 
transmission, especially in weak signal cases. Cyclostationary 
detection is always used to solve weak signal detection, however, 
the computational complexity prevents it from wide usage. In 
this paper, a fast cycle frequency domain feature detection 
algorithm has been proposed, in which only feature frequency 
with significant cyclic signature is considered for a certain 
modulation mode. Simulation results show that the proposed 
algorithm has remarkable performance gain than energy detec- 
tion when supporting real-time detection with low computational 
complexity. 

Index Terms — Cognitive radio, cyclic frequency, cyclostation- 
ary detection, energy detection 

I. Introduction 

The remarkable growth of wireless services over the last 
decade demonstrates the vast and increasing demand for radio 
spectrum. However, the spectrum resource is limited and most 
has been licensed exclusively to users which can work within a 
limited frequency band. Recent study [1] shows that the actual 
licensed spectrum is largely unoccupied most of the time. 
Thus, cognitive radio (CR) has been proposed to solve this 
problem [2], [3]. By sensing and adapting to the environment, 
CR users are able to fill in spectrum holes and serve its 
users without causing harmful interference to the licensed user 
Therefore, the CR system requires spectrum sensing technique 
that detects the unoccupied spectrum band as quickly and 
accurately as possible for its implementation. 

Various spectrum sensing techniques have been presented- 
including matched filter, energy detection and cyclostationary 
detection [4-7]. Matched filter requires perfect knowledge of 
the primary users' signal features. Energy detection method 
is sensitive to noise and interference level. As an alternative, 
cyclostationary detection has a good performance in low 
SNR scenarios, but it costs a large amount of computational 
capacity, which makes it not suitable for real-time detection 
[8]. In this paper, a fast cycle frequency domain feature 
detection algorithm has been proposed, in which only feature 
frequency with significant cyclic signature is considered for 
a certain modulation mode. The detection performance and 
computational complexity for the proposed algorithm are 
analyzed and compared with those of energy detection. 

The rest of this paper is organized as follows: The system 
model under consideration is discussed in Section II. Section 
III gives a brief overview of cyclostationary spectrum analysis. 



Section IV introduces the proposed algorithm. Simulation 
results and performance comparison with energy detection are 
given in Section IV. Finally, conclusions are drawn in Section 
V. 

II. SYSTEM MODEL 

The spectrum sensing problem can be modeled as hypoth- 
esis testing. It is equivalent to distinguishing between the 
following two hypotheses: 



Ho : y{t) = n{t) 

Hi : y{t) = x{t) + n{t) 



(1) 



y{t),x{t)and n{t) denote the received signal, the primary 
user's transmit signal, and the noise, respectively. Hi and 
i?o represent the hypothesis that the primary user is active 
or inactive. Due to the existence of noise, a certain threshold 
Ath should be set to decide whether a primary user is active 
or not. Probability of detection (Pd) and false alarm {Pf) are 
defined to evaluate the detection performance: 



Pd = p{yi > Ath l-ffi ) 
Pf ^p{yi > Ath \Ho) 



(2) 



The goal of detection is to maximize the Pd while maintain a 
given Pf. 

III. CYCLOSTATIONARY SPECTRUM ANAYSIS 

The cyclic autocorrelation of a complex-valued time series 
x{t) is defined by [7]: 
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where a is referred to as the cycle frequency. 

The Spectral Correlation Density (SCD) of x{t), which is 
also known as the cyclic spectrum, can be obtained by Fourier 
transforming the cyclic autocorrelation 

oo 
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If x{t) is Additive White Gaussian Noise (AWGN), then 



when a ^ 0,5"° (/) = 0. 

The frequency-smoothing method is adopted to estimate 
SCD, and the discrete formulation can be expressed as [9] 
[10]: 
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(5) 
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N is the number of input signal samples, x[k]—x [kTs), Tg 



is the time-sampling increment, Fg — Va^T frequency 



sampling increment, I 



//p^J is the discrete value of the 
frequency, and W (v) is the frequency-smoothing function 
centered at u = of width L . 



IV. Cycle Frequency Domain Feature Detection 
Algorithm 

The Spectral Correlation Density (SCD) of a signal is 
a cross-correlation function between frequency components 
separated by / + ^and / — ^, we mapped it from f ^ a 
square to a axis through the following expression 
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And Ix (q) can be regarded as the SCD of a signal at 
the cyclic frequency domain. Different signal modulations 
will exhibit different cyclic features in the cyclic frequency 
domain, as shown in Figure 1 and 2. AM and BPSK 
modulated signals both have significant features at ag — 2fc- 
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Fig. 1. Cycle Frequency Domain for AM signal 

The cyclic features of a signal are decided by the modulation 
types it employed. From the experiment results we found that 
there exists a significant feature at ao = 2/c for both analog 
and digital amplitude modulated signals, where the carrier 




Fig. 2. Cycle Frequency Domain for BPSK signal 



frequency is fc- In the following section of this paper, only 
SCD of the most significant features (ao) of a modulated 
signal will be considered when performing detection. 

We assume that the modulation type and the carrier fre- 
quency fc of the signal are known, namely the cycle frequency 
ao = 2/c is known to us. After SCD estimation at cycle 
frequency ao , the system detection model (1) changes into 
the following form: 



ilo:5,"°(0 = ^r(0 
ili:5^"o(0 = ^r(0 + ^r(0 



(8) 



The cyclostationary spectrum sensing metric of the received 
signal is obtained by 



M = max 
I 



ST (0 



(9) 



A. When the primary user is inactive ( Ha hypothesis), the 
sensing metric becomes 



Afo = max S^" {I) = max \ST (0 I 



(10) 



Theoretically, n{t) is AWGN, so when ao ^ 0, S^J'il) = 0, 
Mo = 0. However, for Limited length SCD, when ao 7^ 0, 
S'^,'{l)^i\ A/o /O [11]. 

B. When the primary user is active {Hi hypothesis), the 
sensing metric becomes 



Ml = max 1 5"" (Z) I 



= max 15^(0 



= max 
I 



max 15^(01 
: max 15^(01 
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(11) 



Mo 



The decision on the presence of a primary user is simplified 
to distinguish Afo from Mi. when the sensing metric of the 
receiving signal is M, namely decide the channel belongs to 



Hq or Hi. a decision threshold Xth is set to classify the metric 
M: 

r Mo(Fo) : M<Xth 

\ Mi{Hi) : if M> Xth ^ ' 

Probability of detection Pd and probability of false alarm 
Pf then be formulated as 

r Pd = Pr(Af> Ath|J?i) 

\ Pf ^PriM > Xth\Ho) ^ ' 

The threshold Xth can be selected for finding an optimum 
balance between Pd and Pf . 

V. SIMULATION RESULTS 

In this section, Monte Carlo simulation results are presented 
to show the better detection performance for the proposed 
algorithm when comparing with energy detection. Simulation 
parameters are listed in TABLE L 



TABLE I 
Simulation parameters list 



Parameter 


Value 


Modulation type 


AM 


Canier frequency 


1 MHz 


Bandwidth 


10 KHz 


Sampling frequency 


3 MHz 


Sampling time 


1.365ms 


Channel 


AWGN 


Window type 


hamming 


Frequency smoothing length(L) 


1300 


Sampled data length(N) 


4096 



Figure 3 shows the curves of receiver operating character- 
istics (ROC) of the proposed method and energy detection, 
it's obvious that the proposed method outperforms the energy 
detection algorithm. 

Energy detection performs badly when SNR bellows -17db, 
while proposed method still has a good performance under 
even -22db. For example, when SNR equals to -22db and 
the probability of false alarm {Pf) is 0.1, the probability of 
detection (Pd) of the proposed method and energy detection 
equal to 0.99 and 0.18, respectively; when Pf = 0.01 , the 
two Pd change to be 0.93 and 0.25, respectively. A cognitive 
radio operation in licensed TV bands (IEEE 802.22 working 
group) defines "required SNR sensitivity" for primary user 
signals to be: -22dB for DTV signals and -10 dB for wireless 
microphones [12]. Therefore, the proposed method meets the 
requirements well. 

If we choose EFT to deal with the formulation (5) and use 
the following expression to calculate the energy of a signal 
(Ex), the complexity of the two methods are shown in table 
2. 




10"^ 10" 10 ' 

Prability of false alann(Pf) 



Fig. 3. Performance comparison between the proposed method and energy 
detection 



N 

(14) 

4=1 

where x{i) is the sample of a signal, N is the sample number 



TABLE II 

Complexity comparison between the proposed method and 
energy detection 





The proposed method 


Energy detection 


Real multiply 


2Nlog2N + 5L 


4N 


Real add 


3N\og2N + 3L 


3N 



From Table 2, we can see that the complexity of the 
proposed method is approximately log2 N times of energy 
detection. The main computation consumption is from EFT 
calculation in the method, for a modern FFT chip this is 
not a big matter as if N is not too big. So the more extra 
computation complexity of the proposed method than energy 
detection is worthy, when considering the performance it 
achieves. To achieve better detection performance, energy 
detection demands much longer sampling time [13], which 
is not necessary for the proposed method. Therefore, the 
proposed Fast Cycle Frequency Domain Feature Detection 
algorithm support real-time detection as well. 

VI. Conclusion 

In this paper, a fast cycle frequency domain feature detec- 
tion algorithm is proposed, in which only feature frequency 
with significant cyclic signature is calculated. Compared to 
the traditional implementation of cyclostationary detection, 
computation complexity for the proposed algorithm is greatly 
reduced under the finite prior knowledge of modulation type 
and carrier frequency. Simulation results show that the de- 
tection performance of the proposed method outperforms that 
of energy detection. Therefore, the proposed method is more 
suitable for real time spectrum sensing in cognitive radio 
system. 
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